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STABILITY ANALYSIS OF MATRIX WIENER HOPF 
FACTORISATION OF DANIELE KHRAPKOV CLASS AND 
RELIABLE APPROXIMATE FACTORISATION 

ANASTASIA V. KISIL 


Abstract. This paper presents new stability results for matrix Wiener-Hopf 
factorisation. The first part of the paper examines conditions for stability 
of Wiener-Hopf factorisation in Daniele-Khrapkov class. The second part of 
the paper concerns the class of matrix functions which can be exactly or ap¬ 
proximately reduced to the factorisation of the Daniele-Khrapkov matrices. 
The results of the paper are demonstrated by numerical examples with partial 
indices {!,—!}, {0,0} and {—!,—!}. 


1. Introduction 

This paper examines the stability of Wiener-Hopf matrix factorisation mmm 
in a certain class of matrices. In the essence, a factorisation of a scalar or matrix 
function G(f) is its decomposition into a product 

G(£) = G+(t)G_(t) (1-1) 

with the invertible factors G+(i) and G_(i) analytically extendable into the up¬ 
per/lower half-plane ('Section l2.Hl . We consider the class of Daniele-Khrapkov 2x2 
matrices BCE], which have the form 

K(t) = I + /(t)J(t), (1.2) 

where f(t) is an arbitrary scalar function with algebraic growth at infinity, and J(t) 
is a polynomial matrix with 

J 2 (f) = A 2 (f)I, 

where A 2 (t) is a polynomial in t, and I is the 2x2 identity matrix. The Daniele- 
Khrapkov matrices can be factorised explicitly (Section []J. 

Fundamentally, the scalar and the matrix Wiener-Hopf factorisations (O) are dif¬ 
ferent: the former has a constructive solution in terms of a Cauchy type integral and 
the latter has no explicit solution in general. The existence of matrix factorisation 
under general assumptions has been proved by Gohberg and Krein |13j . Neverthe¬ 
less, up to date constructive matrix factorisation remains a formidable challenge. 
Due to its complexity, different classes of matrix functions have to be treated sep¬ 
arately (see the recent review article ED- The class of Daniele-Khrapkov matrices 
is very important for applications and arise naturally in a number of interesting 
problems in acoustic, electromagnetics, etc., see for example |T[[gIf2T|. 

The Wiener Hopf factorisation m is said to be stable if small changes in the 
matrix function G (t) lead to small changes in the factors G+(t) and G_(f) (Sec¬ 
tion a. Almost all implementations of Wiener-Hopf technique are performed nu¬ 
merically [l^ , therefore a careful analysis of stability is essential. Among popular 
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approximate techniques are truncated pole removal m and rational approxima¬ 
tions There are also new asymptotic methods HEBEI], which also rely 

on stability. Even in the rare cases when explicit factorisations are known, e.g. 
for Daniele-Khrapkov matrices, they still require numerical computations of scalar 
factorisations. Those computations introduce small errors, which can lead to large 
errors in the Wiener- Hopf factors (Section [3|. 

A landmark theorem of Gohberg and Krein (T5] § 6.2] gives general conditions for 
stability of matrix Wiener- Hopf factorisation (Section[3]). The difficulty of applying 
these results is that the stability conditions depend on the knowledge of Wiener- 
Hopf factorisation and hence are impractical to check. The aim of this paper is to 
provide direct criteria for stability of factorisation in a case of Daniele-Khrapkov 
matrices. The conditions are demonstrated by numerical examples. 

This work is a continuation of the author’s paper [161. which demonstrated a novel 
method of approximately solving scalar Wiener-Hopf equations. In the scalar case 
the formula for the solution in terms of a Cauchy type integral was used to bound 
the error in the factors. In this paper the previous results are extended to the 
Daniele-Khrapkov matrices. 

The first part of the paper establish stability of the Daniele-Khrapkov class under 
perturbations within the class. There are benefits to considering the ‘near’ matrices 
only within the class. It allows to answer the question if numerical implementation 
of the factorisation is stable. This also allowed to obtain explicit error bounds. The 
third advantage is that in a specific case stronger results can be obtained then in 
the general case. 

The second part of the paper extends the class of matrix functions to these which 
can be approximately reduced to Daniele-Khrapkov matrices. The class of matrix 
functions considered by Abrahams in [2] is a special case of this construction. It is 
shown that the stability results could be applied to this meromorphic factorisation. 
This is then used to show stability in a interesting example. 


2. Preliminaries 


Throughout the paper we are using the subscripts + and — to denote functions 
which admit an analytic continuation into the the upper and lower half-planes 
respectively. The Wiener algebra W (R) over the real line [Tj2] Ex. 2.2] consists of 
all complex valued functions / in R that admit a representation of the form 


OO 



— OO 


for some d £ C and k £ L i(R). 


2.1. Wiener—Hopf factorisation. This subsection recalls the different types of 
Wiener Hopf factorisation, which have their own merits, see m for a detailed 
exposition. Let G (t) be in the matrix Wiener algebra W 2 X 2 (R) [HI § 5.2]. If 
det G (t) ^ 0 for all real t then there exists the full factorisation 


G(t) = G+(t)D(t)G_(t), let, 


(2.1) 
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where factors and their inverses belong to the subalgebras of analytically extendable 
functions to the respective half-planes 


G + 1 e W+ 2 ( 


G^ 1 G W 


and D(f) = diag 


t — i 
t + i 


2x2 v 1 

t-F K2 - 


t + i 


The integer exponents k± and «2 are called partial indices. Unlike factorisation, 
the partial indices are unique. But in contrast to the scalar case, they cannot be 
determined a priori in general. 


A factorisation (ED with the invertible factors G+{t) and G_(i) analytically ex¬ 
tendable into the respective half-planes and polynomilally bounded growth at in¬ 
finity will be called function-theoretic factorisation. The function-theoretic factor¬ 
isation is useful in applications since it retains most information and is easier to 
find. 


Remark 2.1. The partial indices are linked to the growth at infinity in function- 
theoretic factorisation, see [Q. 


It is also useful to consider a meromorphic factorisation, where the conditions are 
further relaxed to allow the presence of a finite number of poles and zeroes in the 
factors. 


2.2. Scalar Error Estimates. The index of a continuous non-zero function K(t) 
on the real line is: 

ind(iGi)) = -!-( lim arg K{t) — lim argAf(f)). (2.2) 

Z7T v t—>■+oo t—t—oo 


Note that ind|—| = 1. Thus, given a function K(t) with index n one can reduce it 
to zero index by considering 



For the rest of this subsection it will be assumed that all functions have zero index. 


We also assume that K{t) — > 1 for t —> ±oo, then we can normalise factors such 
that K±{t) —> 1 for t —> ±oo. A non-zero Holder continuous function K(t) on the 
real line with K(t) — 1 in L 2 (R) possesses a factorisation [TO] 

K(f) = K + (t)K_(t), 

where K± ( t ) are limiting values of functions analytic and non-zero in the respective 
half-planes. 


The distinctive feature of the scalar factorisation is the ability to express the factors 
in terms of the Cauchy type integrals. It is the existence of such expressions and 
the bounds in L p on the Hilbert transform which allowed to obtain some useful 
estimation DSI- We adapt them here for L 2 case in the following form. 

Theorem 2.2 (Additive Estimates in L 2 ). Let F(t) = F + (t ) + F_(t) and F(t) = 
F + (t) + F_(t) with ||-F(t) — F(t)|| 2 < e then 

\\F±(t)-F±(t)h<e. 


Theorem 2.3 (Multiplicative Estimates in ^ 2 )- Let K(t) = K + (t)K_(t ) and 
K(t ) = K + (t)K-(t) be two functions and m < |AT| < M. If \\K(t) — K(t) H 2 < e 


then 


\\K ± (t) - K ± (t)h < 


5 (M + e) 1 / 2 

(m ~ e) 6 
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The above results are special cases of theorems from m with some more explicit 
constants calculated. 


3. Stability of Matrix Wiener-Hopf 


For a sake of completeness we review here the most general results on stability of 
matrix factorisation, since they are not widely known in the Wiener Hopf com¬ 
munity. The examples are adapted from a different context of a Riemann-Hilbert 
problem on a circle. There is a wealth of different classes of factorisations con¬ 
sidered by different authors, for the purpose of clear exposition we consider here 
only factorisation in Wiener algebra dm 


The simplest example of instability is obtained by mapping an example m from 
the unit circle to the real line. Consider a diagonal matrix function with partial 
indices {1, —1} 

M 0 \ ( M 0 \ 

o sMf «) L (3l) 

Perturbing the matrix we have 


t—i 

t-\-i 

e 




(3.2) 


This example demonstrates that a small perturbations can not only change the 
factors by an arbitrary amount but can also change the partial indices (from {1, 
— 1} to {0, 0}). This is significant because the partial indices are uniquely defined. 
Note that the sum of the partial indices remains the same. This is true in general, 
which can be demonstrated if we equate the determinants of both sides to reduced 
the problem to scalar factorisation. The partial indices add to give the index (12.211 
of the determinant. In this case, the index of a function / is the winding number 
of the curve (Re /(f), Im /(f)), t € R. Hence, ind(/) and, thus the sum of partial 
indices, are stable under small perturbations. 


Remark 3.1. It is possible to use the non-uniqueness of factorisation m to obtain 
a different factorisation of m 


W 0 

0 B 


1 t—i 

t+i 


Wi ° 

0 w 


1 t + i 
e(t-i) 

0 1/e 


(3.3) 


This is more similar to m- 


The following surprising theorem provides the necessary and sufficient conditions 
for the partial indecies to be invariant under sufficiently small perturbations. 

Theorem 3.2 (Gohberg - Krein, [19] § 6.2]). The system «i > • ■ ■ > K n of partial 
indices is stable if and only if 

Kl ~ Kn < I- 


In fact, this condition is also sufficient for the stability of factors in the Wiener 
norm. 

Theorem 3.3 (Shubin, JT9] § 6.6]). Assume the matrix function G has a Wiener- 
Hopf factorisation and the tuple of its partial indices is stable. Then, for every 
e > 0 there exists a <5 > 0 such that, for ||F — G|| < 6, the matrix function F admits 
a factorisation in which ||F± — G±|| < e. 
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An obstacle in using this result in applications is that one cannot in general de¬ 
termine the partial indices without constructing the factorisation. The next section 
presents new conditions for stability of factorisation for Daniele-Khrapkov matrices. 

4. Error estimates in Daniele-Khrapkov Matrices 

This section examines function-theoretic factorisation of matrices of Daniele-Khrap¬ 
kov class G2D- This class was first considered by Khrapkov in connection to static 
stress fields induced by notches in clastic wedges m ■ There are other numerous 
applications, e.g. related to wave propagation miana. 

Due to this special form (11.211 K (£) can be re-expressed as 

K(i) = r(£) (cosh[A(£)0(£)]I + sinh[A(£)0(£)]J(£)), 

where: 

»<‘> = a^^§§)' (4 - 1) 

Multiplication of the above matrices is commutative, moreover 

Ki(£)K 2 (£) = !?(£)( cosh[A(£)0(£)]I + sinh[A(£)0(£)]J(£)), 

where: 

R(t) = r x (£)r 2 (£), 0(f) = 9i(t) + 9 2 {t). 

This property is enough to obtain function-theoretical factorisation 

K±(£) = r±(t)( cosh[A(t)0±(t)]I + ^ sinh[A(£)0±(£)]J(£)), (4.2) 

where 

r(t) = r-(t)r + (t), 0 = 9-(t) + Q+(t). 

The limitation is the degree of the polynomial A 2 : if it is greater than two then 
cosh[A(£)0±(£)] and sinh[A(£)#±(£)] have exponential growth at infinity [2], This is 
an obstacle to the use of the Wiener-Hopf technique. 

We consider the question of stable factorisation for Daniele-Khrapkov matrices in 
the following sense. Let K (£) and K (£) be of Daniele-Khrapkov type and suppose 
|| K(t) — K(t) || 2 is small. We provide an estimate on || K±(t) — K±{t) || 2 . This 
splits into three parts. The first part is to establish estimates for ||r(t) — fit) || 2 
and ||0(t) — 9(t) || 2 defined by (14.11) . The second is to apply the error estimates 
to parameters r±{t) and 0±{t) of the factors. Lastly \\K±(t) — K±{t) || 2 can be 
examined. 

Consider the matrix function K(t) and its perturbation K (t) 

K{t) = I + f(t)J(t), K it) = I + /(t)J(t), 

such that || A (t)f(t) — A(t)/(t)|| 2 < e. In this setup the perturbation of r(t) can be 
estimated as follows. 

Lemma 4.1. Let r = \J\ — A 2 (t)/ 2 (f) and f = yj 1 — A 2 (t)/ 2 (f). Suppose that 

the winding number of (1 — A 2 (t)f 2 (t)) is zero, then for ||A(£)/(£) —A (£)/(£)||2 <e 
the following estimate holds 

h r l | 2 < —C 

m 

where m = minR{|r(£)|, |f(£)|} > 0 and N = maxR{|A(£)/(£)|, |A(£)/(£)|} < oo. 
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Remark 4.2. The assumptions are natural since |r(f)| 2 is the determinant of the 
matrix K which together with the determinant of its inverse is non-zero. 


Proof. Since winding number of (1 — A 2 (t)f 2 (t)) is zero and e is small enough we 
have winding number of (1 — A 2 (t)/ 2 (t)) is also zero. The square root for r(t) 
in (14.111 can be taken single valued. In the inequality 


| \[a — Vb\ 


\a~ b \ < \a-b\ 

y/a + y/b 2 min( v / a, y/b) 


we substitute a = 1 — A 2 ft) f 2 ft) and 6=1 — A 2 ft) f 2 ft). We also replace 
min(- v /a, y/b) by a smaller value m = minR{|r(f)|, |f(t)|} > 0. Integrating squares 
of the both sides over the real line we obtain 


||r-f|| 2 < i-||A 2 (/ 2 -/ 2 )|| 2 

< -L||(A(/ + /))(A(/-/))|| 2 

< ^ (jf \mm+Kmmm - m \ 2 <**) 1/2 

N 

< —e, 

m 

since \A(t)f(t)\ and \A(t)f(t)\ are bounded by IV = maxR{|A(t)/(f)|, |A(t)/(£)|}. 

□ 


Similarly the behaviour of 6 under perturbation is important. 


Lemma 4.3. Let 


6(t) 


A (t) 


In 


( l-A (£)/(£) \ 

Vi + a wmj' 


m 


A (t) 


In 


( l- A(t)f(t) \ 
\l + A(t)f(t)J ' 


Suppose that the winding number of ( i-|.A(*)/(t) ) ® s zero, then for small || A(t)/(t) — 
A (£)/(£) || 2 < e the following estimate holds 


\\e-~e\\ 2 < 


2e 

cd 2 L ’ 


where d = minR{|l + A(t)/(f)|,|l + A(£)/(£)|} > 0 and L = maxR 
c = min R |A(£)| > 0. 


l—A (t)f(t) 
l+A (t)f(t) 


Remark 4.4. Since A has no zeroes on the real line we can assume min|A| > c > 0. 
Also note that |1 + A(t)f(t)\ and |1 — A(t)f(t)\ are non-zero and finite respectively 
since they are multiples of detA'. 


Proof. From the assumption on zero winding number, the logarithms in the defini¬ 
tion 6(t) and Oft) are single valued functions. The the mean value theorem applied 
to the logarithm function provides an inequality: 


| In a 


In b\ < 


\a-b\ 

min(a, 6) 
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We substitute In a = A (t)9(t), In 6 = A (t)6(t) and replace min(a, b) by L defined 
in the statement. Then, squaring both sides and integrating over the real line we 
obtain: 


\\9-9h < 


< 


< 


1 

cL 


i - a (t)m 

l + A (£)/(£) 


i - mm 

1 + A (£)/(£) 


2 

cL 

2 


A(£)/(£) - A(£)/(£) 


cd 2 L 


(l + A(£)/(£))(l + A(£)/(£)) 

l|A(t)/(i) — A(i)/(t)|| 2 , 


where c and d are defined in the statement. 


□ 


Now we are in the position to apply the scalar error estimates. Under the assump¬ 
tions of the above Lemma 14.31 and using the additive error estimates Theorem 12.21 
we obtain: 

\\e±-~e±h<^ z e. (4.3) 

Using Lemma ITT! and the multiplicative error estimates Theorem l2.3l it follows that 

5 MN 


lr± - r ±|| 2 < 

where M = maxR{|r(t)|, |r(£)|} > 0. 




(4.4) 


To simplify calculation in the next theorem, we will asssume that 


J = 


0 

ko 


ki 

0 


is a constant matrix. Then, a sufficiently small ||A(£)/(£) — A(t)y(t) || 2 guarantees 
that ||K - K|| 2 is small as well. 


Theorem 4.5. Let K and K be of the above form, \\A(t)f(t) — A(f)/(t)|| 2 < e 
and A (t) = C, satisfying the assumptions of Lemmas \ 4.1 \ an d \4-S\ ■ Then, the error 
||K±-K±|| 2 is a linear function of e and the exact estimates can be obtained using 
the above scalar estimates. 


Proof. Let an and an are the top-left elements of K and K respectively. Then 
||on — o-ii|| 2 = ||t±(£) cosh[A(t)0±(i)] — r±(t) cosh[A(£)0±(£)]|| 2 

< ||r±(cosh[A(£)0±(£)] - cosh[A(£)0±(£)])|| 2 

+ || cosh[A(£)0±(£)] (r± - f±)|| 2 , 

where the triangle inequality was used. Then, using the mean value theorem for 
cosh we obtain 

||an — Sn ||2 < |t±| |sinh[A(£)6>±(£)]| ||A(£)0±(£) - A(£)0±(£)|| 2 

+|cosh[A(£)0±(£)]| ||(r±-f±)|| 2 . 

To complete the calculation it is enough to use the bound for |r±|, | sinh[A(£)0±(£)]| 
and | cosh[A(£)0±(£)]|. This follows from r± and 6±, being bounded, having zero 
winding number and tending to a constant |161 . The calculations for other entries 
||aij — dij || 2 , i, j = 1,2 are performed analogously. All the norms of 2 x 2 matrices 
are equivalent so it does not matter which one is chosen. □ 

In the subsequent Sections we present several situations where our results may be 
applied. Numerical example will be presented in Section [6] 
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5. Approximate reducing to extended Daniele-Khrapkov 

5.1. Exact reduction to Daniele-Khrapkov Matrices. The most general class 
of matrix functions which can be factored using the above technique is 

K = S + ( ffl I + g 2 J)S_. (5.1) 

with S + and S_ analytic in the upper and lower half-plane respectively and trJ = 0. 

This can be rearranged as 

K = £iSi+ 52 S 2 , where Si = S+S_ and S 2 = S+JS_. (5.2) 

The challenge is to work backwards from Equation •:E2 to & The first step 
is the factorisation of Si = S + S_ and second step is to ensure the second term 
satisfies the necessary conditions for J = S^ 1 S 2 SI 1 . To satisfy these considerations 
one can take Si and S 2 to be rational, this class was studied in m- 

Now we will outline the procedure to reduce Eauation l5.2l to l5Tl Initially one must 
rule out the case when Si has a zero on the real line. Since the matrix K does 
not have any zeros, any zeros of Si must be compensated either by multiplying by 
fi or by adding / 2 S 2 . So by constructing a different linear combination it can be 
assumed that Si is non-zero on the real line. Then using the rational factorisation 
Si = S + S_ we obtain 

K = S + ( ffl I + ff2 R)S_, 

with R = S + -1 S 2 S_ -1 . 

This can it can be re-written as 

K = Si+(/ 1 I + / 2 J)Si_, 

where J = R — l/2tr(R) for some new functions fi and / 2 , see [23] for further 
details. We will call such matrices extended Daniele-Khrapkov class. 


5.2. Approximate reduction to Daniele-Khrapkov. We give a description of 
a larger class of matrices which may approximately factorised through approxima¬ 
tion by matrix functions from the extended Daniele-Khrapkov class (1521) . Those 
matrices have the property that every entry of the matrix has elements of the form: 

/'i T f'2 e , 

with two fixed arbitrary functions /i and / 2 and rational functions r],- and rfj. 
In the whole generality it shall be discussed elsewhere. Here, we concentrate on a 
subclass, related to work [2] with interesting applications [1. . This subclass allows to 
overcome the problem of exponential growth of the factors in the Daniele-Khrapkov 
matrices for high degree of polynomial A (t). This approximate procedure is simpler 
than the exact one provided by Daniele [El § 4.8.5]. 

Let us begin with matrix 

K<‘> = I + /W (^) t’)' 

We can rearrange it into the form 

K(i)=I + s(t)J(t), 


J(<) 


0 ( 5 @) 1/2 ) 

( S §)‘ /2 o J ’ 


with: 
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and g(t) = fi^ij//) 1 ^- The advantage of this rearrangement being, 

J 2 (f) = I, 

and the disadvantage is that now J has branch cut singularities. To overcome that 
Abrahams proposed to rationally approximate (y^yy) 1 ^ by r/v(i) giving 


3 N it) = 


0 l/r N {t) 
r N {t) 0 

This procedure is exact when n{t) and p{t) have perfect squares as factors. 

The approximate matrix can be decomposed as in (Id.2D 

Kvr(f) = I + g{t)JN{t) = Qn-Qn+, 

but the factors Q./v± have poles. Hence, a meromorphic factorisation is obtained. 


Remark 5.1. Error bounds (TOl) and (14.4D on 9± and r± still hold in this mero¬ 
morphic factorisation. 


To remove poles we can consider the factorisation 

K N (t) = (Qjv_M)(M~ 1 Qjv+), (5.3) 

where M is a rational matrix, which is chosen such that the resulting factorisation 
has no poles in the required half-planes, see [2] for further details. We are turning 
to illustrations of this method. 


Example 1. This example is concerned with the earlier example of instability m- 
The aim is to show that although the indices are 1 and —1 it is still possible to have 
a stable perturbation. The construction is based on the results from the previous 
sections. 


with 


W «/(*) 

«/(*) 


w 0 

o n 


+ e f(t) 


t + i 
0 


0 

i 

t-\-i 


K 


0 1 

c 0 

0 


i 

t—i 

0 t — i 


K = I + e/(t) 


0 (t — i) 1 {t + i) 1 

c(t — i)(t + i) 0 


The matrix K is of Abrahams type with the ratio of the off-diagonal elements being 
a square. Hence, there is no need for rational approximation and the procedure is 
exact in this case. One can construct the factors using the (14.2D . Lemmas 14.11 
and roi can be applied when / satisfies their assumptions. Hence, a meromorphic 
factorisation has been obtained which is stable for small e. Then, the final step is 
to construct a matrix M as in (15.3D . In the case when f{t) = k the matrix M takes 
the form 

, i/2 , i/2k i/2 , i/2k \ 

W^ l+i W l+i j ; (5.4) 


with det M = 1. This completes the factorisation of the perturbed matrix. 


6. Numerical Results 

This section presents two approximate scalar factorisations with different indices 
and these are used to construct two approximate Daniele-Khrapkov factorisations. 
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6.1. Rational approximation. Rational approximation of functions has its uses 
in Wiener-Hopf factorisation. One example was mentioned in previous section. 
Paper [161 applyies rational approximation to simplify the scalar factorisation and 
avoid calculations of a Cauchy type integral. 

Rational approximation is useful for Daniele-Khrapkov factorisation because once 
the approximations for K\ and K 2 are obtained algebraic expressions such as 

Ki + c , K \ + K 2 , A , /y' 2 . 

can be factored easily. This is not true in general as can be seen from the next two 
examples. 

Example 2. Consider the function with zero index 

F ^ = \j t 2 + k 2 ’ ( 6 ' 1 ) 

and with finite branch cuts from i to ki and from —i to —ki. This function is closely 
associated with the matrix function factorisation from problems in acoustics and 
elasticity, see [T. The factors can easily be seen by inspection 

F * (i >= <IWk' F+<t) = F - ( - ty 

However, the factorisation of F(t) + 1 cannot be achieved by inspection. Rational 
approximation of ^j * 2+4 had been also extensively studied in [151. The approxim¬ 
ation was achieved by constructing an appropriate transformation from the whole 
real line to the unit interval. As a result, an approximate factorisation has a small 
global error (10 -12 on the real line). Here, we produce Figure[lJ which demonstrates 
the closedness of approximation on the whole complex plane. 


Example 3. Let us consider rational approximation of the function 


K = 


[t T 2 i)(t T 3 i) 
(t — 2 i)(t — 3i) 


( 6 . 2 ) 


with the index — 1 . Again, the function has been chosen to have the explicit exact 
factorisation 

I (t + 2 i){t + ki) \/{t + 2 i){t + ki) ft — i\ -1 t — i 
Y ki) t + i \t + i) y/{t - 2i){t - ki) 

The function-theoretic factorisation has growth at infinity, making it more difficult 
to approximate. Nevertheless, it can be rationally approximated and the error 
|K — K | is presented in Figure[2] Importantly, the error of the factors \K± — K±\ 
is also small (Figure [3|. For more details on rational approximation of complex 
valued functions see [23] . 


6.2. Numerical Matrix Factorisation. The stability result from Section [4] can 
be used in numerical computations. Two different examples are presented. For each 
example the Daniele-Khrapkov factorisation is computed in two different ways. 
The first method is the direct use of a Cauchy integral to calculate the scalar 
factorisation of r± and splitting 6±. So the initial matrix is exact, the factors 
have errors due to computation of Cauchy integrals. In the second method the 
entries of the matrix are rationally approximated and for this matrix the exact 
Daniele-Khrapkov factorisation is obtained. The matrix is approximate but the 
factorisation of this matrix is exact. The first method will be referred as “exact” 
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Figure 1 . Contour lines for the real and imaginary parts of func¬ 
tion F (EH) and its rational approximation F. They are super¬ 
imposed on a full colour image using a colour scheme developed 
by John Richardson. Red is real, blue is positive imaginary, green 
is negative imaginary, black is small magnitude and white is large 
magnitude. Branch cuts appear as colour discontinuities and co- 
alescent contour lines. Produced using MATLAB package zviz .m. 
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x 10 


-7 


-3 -2 


-1 


x 10 


3 

-7 


Figure 2. Error in approximating function K (16.21) by [8, 8] plot¬ 
ted as real against imaginary part. The accuracy of an approx¬ 
imation is denoted by the size of the disc the curve is contained 
in. 



Figure 3. Error of factor K± on the real line plotted as real 
against imaginary part. The accuracy of an approximation is de¬ 
noted by the size of the disc the curve is contained in. 


and the second one as “approximate” although the reader should note that both are 
approximate factorisations. The results of these two methods are then compared 
for each example. 

The first example is 


t 2 + 1/0 1 \ 

i 2 + 4 V -2 0 ) ' 


K!(i)=I + 









STABILITY ANALYSIS OF MATRIX W-H FACTORISATION OF DANIELE-KHRAPKOV 13 


The ideas is to rationally approximate by /at- Then the factorisation of: 


( 1 fN 

V cf N 1 


is computed and compared with the “exact” factorisation. The advantage of such 
an approximation is that there is no need to use the Cauchy formula to find r± 
and 0±. Note that the approximate matrix has all rational entries and hence in 
theory factorisation can be achieved using methods for rational matrix functions. 
But in practice the implemented procedures are unstable, making it impossible. At 
present, very few implemented algorithms Wiener-Hopf exist. For example, there 
has been some attempts recently [6] to produce numerical factorisation algorithms 
for rational matrix functions and numerical algorithms for Riemann-Hilbert prob¬ 
lems mnm- 

The second example is 


K 2 (i) = I- 


(t H- 2i)(t + z) / 0 1 


V -2 0 


Similarly the approximate factorisation is considered by approximating \J 


(t+2i)(t+i) 




Figure 4. The modulus of K\ + and K 2 + on the real line. 


The difference in behaviour on the real line of the two examples can be seen in 
Figure [I] This is because their partial indices are different. The first example 
have partial indices {0, 0} and the second {—1, —1}- These partial indices can be 
computed using the following identity 


1 

cf 


f 

1 


1 1 \ f 1 + c 1 / 2 / 0 

c l/2 _ c l/2 J [ 0 1 - C 1 / 2 / 


1 

2 

1 

2 



The errors are compared in Figure [5] and Figure [6] It should be noted that the 
calculation of “approximate” factors took significantly less computational time than 
the “exact” factors. Besides the natural difference in magnitude of errors (due to 
the difference in errors of rational approximations) the shape of the curves are 
dramatically different. It seems the error in Figure [5] is random and in Figure [G] is 
systemic. This suggests that in the first example the error in “exact” factorisation is 
greater than “approximate” factorisation. So the accumulated errors in computing 
Cauchy integrals is greater than the error in once approximating entries of the 
matrix function. The reverse is true in the second example. 
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Figure 5. The modulus of the difference in an elements of “ex¬ 
act” and “approximate” factors for K t . 



Figure 6 . The modulus of the difference in an elements of “ex¬ 
act” and “approximate” factors for K 2 . 
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